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CN '. Abstract 

^f^ I We analyze the contributions to the neutrinoless double /? decay (Oi^/3/5- 

•/^ ■ decay) coming from the Grand Unified Theory (GUT) constrained Minimal 



Supersymmetric Standard Model (MSSM) with trilinear R-parity breaking. 



^^ . We discuss the importance of two-nucleon and pion-exchange realizations of 



the quark-level Oiy(3(3-decay transitions. In this context, the questions of reli- 
ability of the calculated relevant nuclear matrix elements within the Renor- 



o ■ malized Quasiparticle Random Phase Approximation (pn-RQRPA) for sev- 

jiv ■ eral medium and heavy open-shell nuclei are addressed. The importance of 

D . gluino and neutralino contributions to Oi^/3/3-decay is also analyzed. We review 

'"^ ' the present experiments and deduce limits on the trilinear R-parity breaking 

parameter X'm from the non-observability of Oz^/?/3-decay for different GUT 
/\ • constrained SUSY scenarios. In addition, a detailed study of limits on the 

j^ ■ MSSM parameter space coming from the B -^ Xsj processes by using the 

recent CLEO and OPAL results is performed. Some studies in respect to the 
future Oz^/3/3-decay project GENIUS are also presented. 
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I. INTRODUCTION 

The neutrinoless double beta decay (Oz//5/5-decay) is forbidden in the Standard Model 
(SM) since it violates lepton number by two units (AL = 2). Therefore this decay is a 
sensitive probe for different aspects of physics beyond SM. (For recent reviews see e.g. |I|J^). 
Many generalization of the SM admit the violation of laws of the SM to a small extent. In 
this context, non-observability of the Oi//3/3-decay is used to constrain different extensions 
of the SM like those with the left-right symmetry [|],|[, leptoquarks 0, R-parity violating 
supersymmetric (-^pSUSY) models [|6HT^ and composite neutrinos [|T7i|T8 |. 



The most widely discussed case in literature has been the upper bound on the light 
effective electron neutrino Majorana mass {mj) deduced from the experimental lower limit 
on the half-life of 0i//9/9-decay []T]J^. Currently, the most restrictive limit on {mj) is found 
from Oi//3/3-decay in '''^Ge by the Heidelberg-Moscow collaboration [|l^: {mj) < 0.4 — 1.3 
eV [|I|. The uncertainty of this parameter is due to the ambiguity of 0i//?/3-decay nuclear 
matrix elements. It is expected that the future double beta decay experiment GENIUS pi 



based on 1 tonn of enriched ^^Ge would reach the sensitivity for {rrii^) (0.01 - 0.001 eV). 

Besides the simplest and the best known mechanism of lepton number violation based 
on the mixing of massive Majorana neutrinos advocated by different variants of the Grand 
Unified Theories (GUT) the R-parity violation proposed in the context of minimal super- 
symmetry extensions of the SM (MSSM) is becoming the most popular scenario for lepton 
number violation (see e.g. reviews [0,^). The R-parity is a multiplicative quantum num- 
ber defined as R = (_ i)3-B+l+2S ^j^]-^ g^ l and S being the baryon, the lepton and the 
spin quantum numbers of the particle, respectively. Thus, all the SM particles have R=-|-l, 
while their superpartners have R=-l. We note that the R-parity conservation, which guar- 
antee the baryon and lepton number conservation, is not required by gauge invariance or 



supersymmetry and might be broken explicitly or spontaneously at the Planck scale |16 . 

The R-parity can be broken by involving the bilinear and trilinear terms in superpo- 
tential of the MSSM. The bilinear terms generate non-zero vacuum expectation value for 
the sneutrino field, which leads to the lepton number violating interactions based on the 



neutralino-neutrino and chargino-electron mixing |jl2[. The trilinear terms represent the 
interactions, which violate directly the lepton number and lepton flavor ||^12|. Both ways 
influence the low-energy phenomenology and therefore using some exotic processes like the 
Ou (3 jS-decaj one can impose limits on the parameters connected with the new physics. 
Supersymmetric models with R-parity non-conservation (^jjVISSM) have been extensively 



discussed in the literature (see e.g. pT| , |22[| ) and were also used for the study of R-parity vi- 
olating trilinear term contribution to the Oz//?/5-decay by Mohapatra [§] and Vergados 0. 
The first calculations were concentrated only on the conventional two-nucleon mode of Oi//3/3- 
decay, which assumes direct interactions between quarks of two-decaying neutrons H^-^] . A 
detail study of this mechanism was carried out in Ref. [§]. By using viable phenomenological 
assumptions about some of the fundamental parameters of the ^pMSSM (e.g. the ansatz 
of universal sparticle masses and that the lightest neutralino is bino-like) it was found that 
the limit on the R-parity violating first generation Yukawa coupling A'^^^^ derived from the 
observed absence of the Oz//5/5-decay is more stringent than the corresponding limit expected 



from the forthcoming accelerator experiment at HERA. The authors end up with the conclu- 
sion that the gluino exchange R-parity violating mechanisms of the Oz//5/3-decay dominates 
over the neutralino ones |^. 

Another scenario for reduction of a number of supersymmetric parameters associated 
with the limit on A'^^^ has been outlined in Ref. |p. The authors implemented relations 
among the weak scale values of all parameters entering the superpotential and the soft 
SUSY breaking Lagrangian and their values at the GUT scale. The obtained results showed 
the importance of the neutralino exchange mechanism. Similar studies have been performed 

IM. We note that the GUT's con- 



also within gauge mediated SUSY breaking model 

strained SUSY scenarios have advantages of fewer degrees of freedom and usually predict 

more stringent limits on X[i^ 



23|J24 



Recently, the dominance of the pion-exchange R-parity violating Oz//?/?-decay mechanism 
(based on the double-pion exchange between the decaying neutrons) over the conventional 
two-nucleon one was proven for ^^Ge isotope |T^. We note that the attention to the pion- 



exchange Oz//3/3-decay mechanism was paid out first by Pontecorvo [^ and that this mech- 
anism has been found less important in the presence of light Majorana neutrinos and more 
important, if heavy Majorana neutrinos are considered [^. There are two advantages, 
which favor pion-exchange R-parity violating mode over the two-nucleon mode. First, the 
effective radius of the two-nucleon R-parity violating interaction is small because of exchange 
of heavy SUSY particles and therefore this mode is suppressed by the nucleon-nucleon re- 
pulsion at short distances. On the other side the pion-exchange mode leads to a long-range 
nuclear interaction, which is significantly less sensitive to short-hand correlations effects, 
as is mediated by light particles. Secondly, the enhancement of the pion-exchange mode 
has an origin in the bosonization of the n~— > vr"*" -|- 2e~ vertex and is associated with 
the pseudoscalar hadronic current structure of the effective R-parity violating Oz//?/5-decay 



Lagrangian on quark level ||T^,[T2|. The new pion-exchange mechanism of Oz//3/3-decay offers 
more stringent limits on A'^^^^^ |T^. The one- and two-pion exchange realization of this mode 
have been discussed for experimentally interesting nuclear systems in Ref. [O . Both gluino 



and neutralino exchange mechanisms have been found of comparable importance within the 
phenomenological scenario of this mode. 

The procedure leading to the final constraint on X[^^ consists of two parts: In the first 
step, the relevant SUSY parameters are determined within a proper SUSY scenario. Next, 
the reliable evaluation of Oz//3/3-decay nuclear matrix elements has to be performed. Since 
the importance of the particle-particle residual interaction of the nuclear Hamiltonian for the 
description of open shell nuclear system was discovered [0, the proton-neutron Quasipar- 
ticle Random Phase Approximation (pn-QRPA) has been widely chosen in the calculations 
of the nuclear double beta decay transitions [p],p|,|28|-[35[] . However, the extreme sensitivity of 
the calculated nuclear matrix elements as well as the collapse of the QRPA solution in the 
physically acceptable region of particle-particle strength address some questions about the 
predictive power of the obtained results |]36| , |37| . The renormalized pn-QRPA (pn-RQRPA) 
without [BH] and with proton-neutron pairing E^ 



which take into account the Pauli exclu- 
sion principle, do not collapse and offer more reliable results in respect to the QRPA |]36| , |37| . 
This method has been used also in the recent nuclear structure studies of the Oz//?/3-decay 
matrix elements [p],p!0|- [T^j36[| . Some other extensions of RQRPA can be also addressed in 
this context [1131. 



The goal of the present paper is to perform a comprehensive analysis of the R-parity 
violating trilinear term contribution to 0z//9/3-decay within the GUT constrained MSSM. We 
shall focus our attention on a detailed pn-RQRPA study of the relative importance of the 
pion-exchange and two-nucleon modes of this process for experimentally interesting Qvj3j3- 
decay nuclei. Our aim is also to discuss the present as well as possible future (imposed by 
GENIUS experiment) limit on the trilinear breaking SUSY parameter A'^^^^ and identify this 
limit with the gluino or the neutralino exchange mechanisms. Moreover, all the calculations 
are performed within MSSM constrained by different processes, among which the most 
interesting are limits coming from B — > X57 processes obtained by using the recent CLEO 
and OPAL results. 

The paper is arranged as follows. In the next section (section II) we introduce the basic 
elements of the MSSM with the explicit R-parity breaking. Here we also discuss the main 
formulae relevant for the trilinear R-parity breaking contribution to 0z//9/9-decay and the 
realistic pn-RQRPA nuclear structure method, which will be used for evaluation of nuclear 
matrix elements of interest. In section III we calculate the Oz//5/9-decay matrix elements 
within the pn-RQRPA for experimentally interesting isotopes and analyze their uncertainties 
in respect to details of the nuclear model. In addition, we discuss the parameters of the 
MSSM and the importance of the gluino and the neutralino Oz//5/5-decay mechanisms. The 
experimental constraints on the Oz//5/5-decay are then used to constrain the first generation 
of lepton number violating Yukawa coupling constant Am of supersymmetric particles. We 
close with a short summary and conclusions (section V). 

II. THEORY 

A. R-parity violating minimal supersymmetric standard model 

The MSSM is based on the same gauge group as the SM, It is particle content required 
to implement supersymmetry in a consistent way is the minimal one . It is described by the 
superpotential which in the case of R-parity non-conservation contains both the R-parity 
conserving and breaking parts 

W = W^ + Wft, (1) 

where 

Wo = hlQiH^u] + hf^Q.Hj] + h^jUHde] + fiHdHu. (2) 

and 

Wfi = KjkLiLjul + Xijj^LiQ jdl + \,ijjji1<fjd1 + fijLjHu- (3) 

Here Q, L denote the quark and lepton SU(2) doublet superfields, u^, d'^, e^ the corresponding 
SU(2) singlets and Hu-, Hd the Higgs superfields whose scalar components give mass to up- 
and down-type quarks and leptons. In the R-parity breaking part we set \ijk = \'ljk = to 
avoid the unsuppressed proton decay. 

Since supersymmetry in the low-energy world is broken, the Lagrangian of the theory is 
supplemented with the "soft" supersymmetry breaking terms: They are given by: 



and 



-Csoft = [Af^.^h^Q.H^u'; + Af.^^hf^Q^Hd'j + Af^^^hf^LiH^^ + H.c. 
+ Bi2 {HdHu + H.c.) + mjfjHdl^ + mljH^\^ 

+ {^M,i,Bi^B + ^M2^^^^ + ^m^^^>^" + H.c.) (4) 



-^soft^ — ^ijk^i^j'Ui. + A^ji^LiQ jd^ + A^ji^U^djdf^ + ^2j^j^u + /^Ij-^jr-nrf. (5j 



Here, the fields with tilde denote the scalar partners of the quark and lepton fields, while 
the ipi are the spin-| partners of gauge bosons. 

In the present paper we concentrate on the trilinear terms R-parity breaking superpo- 
tential only, leaving complete treatment of bilinear terms for future studies. 

Since our main goal is the -^pMSSM description of the 0i//9/9-decay, we have to derive 
the effective Lagrangian for that process from the superpotential. We will show below that 
this Lagrangian depends on many free supersymmetric parameters, like masses or coupling 
constants. In order to reduce the number of these parameters, we apply conditions motivated 
by supergravity theories unification for soft mass parameters at the GUT scale (see e.g. 
pO|-^). We shortly present our method of finding the low-energy spectrum of the MSSM 
below. 

We begin our procedure with running the gauge and Yukawa coupling constants up to 
the point where they unify (GUT scale). The initial conditions for the gauge couplings are 
set up by the values of aem(^z) = 1/127.89 and as = 0.118. The values of Yukawa couplings 
at mz are connected with the standard model quark and lepton mass matrices in the usual 
way: 



f sin/? 

71" 



Mu = "-^S''-Y^S''-\ 



72 ^' 



Mo = ^^^^S^-Y^S' 



Me = ^^^S^-Y^S^-l (6) 

v2 

Here, the unitary matrices S are connected with mixing among the MSSM fields of the 
matter sector induced by the electroweak symmetry breaking and K = S^^S^^'^ is the 
Kobayashi-Maskawa matrix. In the matter sector primed mass eigenstates are related to 
their unprimed gauge eigenstate counterparts. 

u' = S^^u + S^^CIF^, 
d! = S^^d- S'^'^C'df, 
e' = S^^e - S^«Ci^^, 
u' = S^^u. (7) 

In above formulas, the following convention for the MSSM up an down-type Higgs vacuum 
expectation values are used: 



(^.°) = ^, (H'u) = ^. v = ^vl + vl = 2AQGeV. (8) 

During the running procedure, the renormahzation group equations (RGE's) change 
between the weak and GUT scales at each particles mass threshold due to decoupling of 
states at the scales below their masses. At first running, we initially set this threshold at 
MsusY = ITeV, using the SM 2-loop RG equations below, and the MSSM RGE's above that 
scale. When the gauge and Yukawa couplings get the GUT scale, we set all the soft scalar 
masses equal to the common scalar mass ttiq, the soft trilinear couplings A to common Aq 
and take gaugino masses to mi/2. In the next step one runs everything down back to the mz 
scale. We want to stress at this point that we do not use the full set of RGE's appropriate 
for the R-parity broken MSSM |H3[ . We estimate that an influence of the R-parity breaking 



constants on other quantities running is marginal due to the smallness of A's. So, we limit 
our attention to RGE's for MSSM with R-parity conserved p^). 



It is well known that running of Tn'j^^ is dominated by negative contribution from the 
top Yukawa coupling, which drives this parameter to a negative value at some scale causing 
dynamical electroweak symmetry breaking (EWSB). Thus the RGE's improved supersym- 
metric potential naturally breaks SU{2) x f/(l)y to f/(l)cm; which additionally allows to 
express some GUT-scale free parameters in terms of low-energy ones. The tree-level Higgs 
potential has the form 

2 2 

Vo-m^\h^\ +m^\hj + m^{h^h^ + ti.c.}^ {\H^\ - \lij ) , (9j 

where ml 2 = tn'jj +IJ''^, rn^ = B/j,, and the phases of the fields are chosen such that m| < 0. 
Due to Q scale dependence of the soft Higgs parameters Vq also strongly depends on Q. Thus 
minimization of Vq can produce very different vacuum expectation values Vu, Vd depending 
on Q. That's why it is necessary to minimize the full one-loop Higgs effective potential. The 
procedure leads to a set of two equations. 



m| ^ jmjj^ + Srf) - {ml^ + S„) tan^ (3 
2 tan^ (3 - 1 






Su, Srf are given e.g. in Ref. [|^. In order to minimize the stop contribution to the finite 
corrections, we take the minimization scale Qmin equal to the square root of mean square of 
stop masses. 

The EWSB causes mixing among many particles. In particular, four gauginos mix to pro- 
duce the so-called neutralinos with the mass matrix M^, which in the ip = {B'^, W^, H^, H2) 
basis has the form: 



(12) 



with Cyj = cos{9w),Sw = sm{9w),S/3 = sin(/3),C/3 = cos(/3) and Mi,M2 being U(l) and 
SU(2) gaugino masses. At this point it is worth mentioning that the physical gluino mass is 
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-nizspcw 
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related to the renormalized DR running mass m^. For very heavy quarks this dependence 
can read PB 



m 



pole 



ms 



l + ^(l5 + 12/( 



(13) 



where m^, and a^ are taken at m^, and the loop function I{r) = ^Inr + |(r — l)^ln(l 
1 for r > 1 with r = rn^/ml. 
Mixing in the gaugino sector effects with four physical neutralinos Xi- 



r ^) + fr 



Xt 



T.^^A. 



V 



1,2,3,4). 



(14) 



i=i 



The matrix A/", diagonalizing the matrix M^ is real and orthogonal. Thus, with real Mij, 
neutralino masses are either positive or negative. If necessary, a negative mass can always 
be made positive by a redefinition of the relevant mixing coefficients Mij — > iMij. 

Similar mixing appears in the slepton and squark sector. The SUSY analog of mixing 
(0) in the SM sector is described as follows: 

-S^^d''* ) ' 



u' 



d' 



e 



^D 



V 



YN^El 



y. 



(15) 



The 6x6 squared mass matrices for the squarks and sleptons have a much more compli- 
cated form and involve parameters from both the supersymmetry breaking and conserving 
Lagrangians. They can be written in the following form: 



Mi 



(3-4sin2^VK)cos2/? 



I '"I 

""" 6 



-/i*M(7COt/? 

V sin [3 gU f, ^^yT (^U J \ 



V2 



-fiMu cot /3 \ 

V sin f3 oUlAU* oUr-\ 
V2 



H — ^ sin^ 9w cos 2/3 



pi/t 



(16) 



/ 5^^m2,5^^t + M2 



Mj 

d 



^D 



D 



-^(3-2sin2^vi^)cos2/3 



—/iMo tan /3 \ 

igDLji^D*gDRJ 



^/l 



Mi 



-Ij,*Md tan (3 

_vcosj3gDRj^D^gDL-\ 

-^(l-2sin2^H/)cos2/3 



S^^ml' 5^«t + M2 



D 



- H* Me tan p 

vcos^^Ejij^E'^ I^E^-l 



■ sin 9]y cos 2/5 



—fiME tan P \ 

vcosfi cEl aE* cEb.] 

^^«m?;' S^«t + Ml 



-m 



I sin^ 9w cos 2/3 j 



pDt 



p£;t_ 



(17) 



(18) 



When the physical masses of scalar particles and tree level values of /i and Bfi are found, 
one can set up thresholds for RGE's of the gauge and Yukawa couplings and repeat the 
procedure by running everything up to the GUT scale. Setting the unification conditions 
for soft parameters again, it is possible to run everything back to the mz scale. At that 
point we repeat calculation of mass eigenstates, run everything up to Qmm and minimize 
the one-loop corrected Higgs potential mentioned above. Repeating such a procedure leads 
to fixed values of n and Bfi. At the end, we check for conditions necessary for existence of 
a stable scalar potential minimum: 



2B^i<2\^i\^ + ml^ + ml^. (19) 



B. Constraints 

Having obtained the low energy spectrum of the model, we are at a good point to test 
our scenario using constraints due to rare Flavor Changing Neutral Currents processes. It 
is well known that FCNC processes may serve as a strong constraint on the supersymmetric 
scenarios. Strong experimental suppression of FCNC transition puts upper bounds on vari- 
ous entries of the sfermion mass matrices at low energy. Many analyses of such constraints 
were performed for the SUGRA MSSM [^. Some of them are also available for the Gauge 
Mediated Supersymmetry Breaking (GMSB) |^. Without going into details of calculations, 
we summarize the basic ingredients of the applied procedure. 

In our considerations we concentrate on limits on the low-energy spectrum coming from 
the B -^ Xs'~f decay. This process is described by an effective Hamiltonian of the form 

mm-- 



-neff = 



-^KlK,,Y.C^{^^)P^{^^), (20) 

where K is the Kobayashi-Maskawa matrix. Pi are the relevant operators, and Ci{fi) are 
their Wilson coefficients. The coefficients C7 and Cg relevant in the subsequent analysis get 
contributions from both the SM and the MSSM interactions. The leading order and next- 



to-leading order SM contributions are discussed in [§^. For the MSSM case, only leading 
order contributions to those coefficients are available |^0 Q. 

The constraints on the low-energy spectrum of our model are connected with present 
constraints on the R-/ parameter which measures the extra (MSSM) contributions to the 
B — > Xs7 decay. Its definition reads: 






^7 = 1+ ^(o)5M. T' (21) 



^NLO QCD corrections for SUSY scenario with charginos and one of the stops lighter than other 



sparticles can be found in |51]. 



where index (0) stands for the leading order (LO) Wilson coefficients, and superscript extra 
for the SUSY (namely charged Higgs, chargino, neutralino and gluino) contributions. 

Limits on the allowed values of Rj are extracted from the present experimental lim- 
its on the branching ratio BR(i? -^ ^sl) recently measured by CLEO collaboration ||51|| : 
BR{B -> X,7) = (2.50 ± 0A7stat ± 0.39sj,st) x 10^^ The same process measured at ALEPH 
collaboration results with BR{B -^ X^-f) = (3.11 ±0.80,tat±0.72,j^,t) x 10"^ p|. Using the 
theoretical dependence of the theoretical expectations of BR(i? -^ -^s7) on Rt, comparison 



with the experimental data |^0[ and taking into account both theoretical and experimental 
errors, we end up with the following estimation of the i?7 range : 



- 6.6 < i?7 < -4.4 or 0.0 < i?7 < 1.3. 



(22) 



In order to find numerical values of R^ given in (pT]) , we use the expressions for C^ 



and Cj 
region 



{0)SM 



listed in e.g. ||4^. Then we exclude points which do not lie in the allowed 

Xg'-f constraints, the free parameter space is limited by 



I). In addition to B 

the conditions of proper electroweak symmetry breaking (|T9|) and the condition of positive 
(mass)^ of mass eigenstates. In the last section we present the resulting constraints in 
tan/5 — rriQ and tan/5 — mi/2 planes for fi of different signs. Moreover, the detailed analysis 
of i? — i> Xsj constraints is shown. 



C. R-parity violating neutrinoless double beta decay 

In this Subsection the main formulae relevant for R-parity violating Oi'f3f3-decaj are 
presented. A detailed derivation of the effective ^p quark-level Lagrangian for the Ou^f]- 
decay and its hadronization additionaly to derivation of the corresponding half-life time of 
this process can be found elsewhere Refs. JSHTTI. 

The trilinear R-parity violating lepton part of the interaction Lagrangian of the MSSM 
takes the form |H|: 



Cx' 



-'^iii 



-Kll[i^L,dL) 



""^ ^d% + {e-L,i7L)dnl^^A + 



-u 



R 



iu-L,dL)dn{ "^^ +H.C.]. 



(23) 



Together with the R-parity conserving MSSM Lagrangian describing interactions among 
gluinos, neutralinos, fermions and sfermions |S^ after integration out of heavy degrees of 
freedom and carrying a Fierz transformation, one can obtain finally the ^pSUSY induced 
quark-lepton interaction for Oz//?/3-decay. 



f-ALe 



2m„ 



e(l + 75)e^ 



'^PS J PS J PS — -VT Jt^Jt^iu 



4 



(24) 



The color singlet hadronic currents in Eq. ( |2^ ) are Jps = u'^'^r>da + u°'da, J^ = u°'a^'^{l + 
l^)dai where a is a color index and a^^ = {i/2) [7^^, 7^^]. The effective lepton-number violating 
parameters rjps and tjt in Eq. (|2^) accumulate fundamental parameters of the ^pMSSM: 



r]ps = Vxe + %f + % + ^9 + 7r/^, 



(25) 
(26) 



with 
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I L 



^— V TTip 



. -I my. 



[^Ri{d)eLi{e) + 
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SL 



eLi{u)tLi{e) 



771 



dR 



(27) 



In Eq. (p^) G^;^ is the Fermi constant, m^ is the proton mass, 0^2 = ^fg/l^Ti^) ^-^d a^ = 5'|/(47r 
are SU(2)l and SU(3)c gauge couphng constants respectively, ms^ , m 



"Mi; 



m.g and ?7i^. 



are 



masses of the u-squark, d-squark, gluino and neutrahnos. We note that the matrix Mij (see 
Eq. (|l^) rotates the 4x4 neutrahno mass matrix in Eq. (|l^) to obtain the diagonal form 
Diag[m,y\. We used the neutralino couplings in the form proposed in Ref. ||53|| : 

elM = -T^{<P)M2 + tanew[T^{(P) - Q (0)1^1, 

eRX4>) = QicfytanOwMii. (28) 

The effective ALg = 2 Lagrangian ( p4D contains contributions from both gluino (SUSY 
parameters rjg and rj'g) and neutralino (SUSY parameters r/^, rj ? and rj^e) exchanges. The 
relevant Feynman diagrams associated with gluino g and neutralino x contributions to the 
Oi//5/3-decay have been discussed in Refs. [|l],||. 

The ^pMSSM model gives the underlying transition of a down-quarks to an up-quarks 
{dd ^ UU + 2e~) only, and results in transformation of neutron into a proton. Till now three 
possibilities of hadronization have been considered. The most natural way is to incorporate 
the quarks in the nucleons which is the well-known two-nucleon mode [§-§|. But the inter- 
mediate SUSY partners are very heavy particles. Therefore in the two-nucleon mode the two 
decaying neutrons must come very close to each other, which is suppressed by the nucleon 
repulsion. Another possibility is to incorporate quarks undergoing the -^pSUSY transition 
not in nucleons but in virtual pions P| JT(]|JT^ . If only one of the initial quark is placed in an 
intermediate pion, we end up with the one pion-exchange mode |T^. If all quarks are placed 
in two intermediate pions, one obtains the two-pion exchange mode, which dominates for the 
Ou (3 jS-decaj of ^Ge ||TO|Jl^ . The possibilities of incorporating quarks in heavier mesons have 
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not been considered so far. However, this way of hadronization is expected to be suppressed 
due to heavier masses of exchange particles. 

The half-hfe for the neutrinoless double beta decay regarding the above three possibilities 
of hadronization of the quarks can be written in the form [jl],|l0|,|l2 



[T,%(0- 



0+)]-^ 



G 



01 






ttN 



where 
M 
M 



2N 
9 

2N 
f 



CA 

ca 

CA 



a 



(0) 



(0) 
V 
V-f 



MfN + dA-qMoTN 
(0) 



vlqMp' + a^A-gMcT' + ar-qMr' 



r4 



a^" (McT-in + Mt-1^) + a'^ (Mgt-2. + Mt-2^ 



(29) 

(30) 
(31) 

(32) 



with c^ = rn?^/{mpme). Here Gqi is the standard phase space factor (see Ref. Il3|j35||) and 
niA = 850 MeV is the nucleon form factor cut-off (for all nucleon form factors the dipole 
shape with the same cut-off is considered). The nucleon structure coefficients a's entering the 
nuclear matrix elements of the two-nucleon mode calculated within the non-relativistic quark 
model (NR) and the bag model (BM) are given in Table 3 of Ref. [|l| . The structure coefficient 
of the one-pion a^'^ and two-pion mode a^'^ are ITDIO: a^'^ = —0.044 and a'^'^ = 0.20. The 



partial nuclear matrix elements of the ^p SUSY mechanism for the Oi^/5/3-decay appearing 
in Eqs. (|Sg)-(|2|) are: 

i+3 



«J 



M 



(0;i E r+r; ^ F,(x) a, ■ a„ |0+), 



«7^J 






Mk = (0;| E rtr; -^ Fk{x) S., |0+), 



(33) 



^3 



where / = FN, F\ J = GTN, GT', GT-kn and K = T', T-kn {k = 1, 2). The structure 
functions Fj^j^k can be expressed as: 



Ffn{x) 
Ff'{x) 

Ft'(x) 



Fgtn{x) 



X 



^(3 + 3x^ +<)e-^^. 



X, 



48 ' 



x2)e ^A, 



Fgt-\A^^) = e ''\ 
Ft-\-k{xt,) = (3 + 3x^ + x\)- 



X: 



2 ' 



FGT-2TT{xTr) 
FT-2nM - 



-■ {x^ - 2) e" 
ix^ + 1) e-^ 



(34) 
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and we used the notations: fjj = (ri — rj)/|ri — rj|, rij = |ri — rj|, x^ = rn^rij, x-,, = niT^rij 
and Sij = 3(jj • Vij aj ■ Vij — ai ■ aj. Here ri is the coordinate of the i-th nucleon and R = tqA^^^ 
stands for the mean nuclear radius with rg = 1.1 fni. 

In order to deduce constraints on the lepton number violating parameters rjps and tjt 
from the non-observability of Oz//5/5-decay it is necessary to evaluate the nuclear matrix ele- 
ments of two-nucleon (A^?^, A^^^) and pion-exchange {Ai"^^) modes within an appropriate 
nuclear model. 

III. NUCLEAR MODEL 

We calculate the nuclear matrix elements within the proton-neutron renormalized Quasi- 
particle Random Phase Approximation (pn-RQRPA) ||38|j39|l , which is an extension of the 
pn-QRPA [0,^] by incorporating the Pauli exclusion principle for the fermion pairs. The 
advantage of the pn-RQRPA over the usual pn-QRPA is no collapsing RQRPA solution 
as one increases the strength of the particle-particle interaction within its physical values. 
Thus the results obtained by the pn-RQRPA method are more reliable, but the pn-RQRPA 
method requires coupled non-linear equation solutions, instead of the usual eigenvalue prob- 
lem in the pn-QRPA formalism. 

The pn-RQRPA method is suitable to deal with the nuclear structure aspects of beta 
transitions of open shell systems and allows to perform calculations in realistic large model 
spaces, which are unaccessible for the shell model calculations. 

The pn-RQRPA formalism consists of two main steps: (i) The Bogoliubov transformation 
smears out the nuclear Fermi surface over a relatively large number of orbitals and (ii) the 
equation of motion in the quasiparticle basis determines then the excited states. 

If only proton-proton and neutron-neutron pairing is considered the particle (c:^^^ and 
c^jn^, T = p,n) and quasiparticle (a;^^^ and a^^^, r = p,n) creation and annihilation 
operators for the spherical shell model states are related to each other by the Bogoliubov- 
Valatin transformation: 



Ur —V-r 
Vr Ur 



(35) 



where the tilde indicates the time-reversal operation (irmr = ("~l)''^~'^^ctT-mr- "^^^ occupa- 
tion amplitudes u and v and the single quasiparticle energies Ej. are obtained by solving the 
BCS equation. Then one gets a nuclear Hamiltonian in quasiparticle representation 

H = ^ Era^^^arm-r + H22 + H^q + Hq^ + H^i + if 13, (36) 

where Hij is the normal ordered part of the residual interaction with i creation and j 
annihilation operators (see e.g. Ref. [0). 

In the framework of the pn-RQRPA the im}^ excited state with the angular momentum 
J and the projection M is created by a phonon-operator Qjm^^ 

|m, JM) = QtiA^lpj,) with QTm^ |0+p^) = 0. (37) 



Here |0^p^) is the ground state of the initial or the final nucleus and the phonon-operator 
Q^M-K is defined by the ansatz: 
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Q7m^ = E^(;n,J^)^HP^, JM) + l^(;„,j.)i(pn, JM). 



(38) 



pn 



A'^{pn, JM) (v4(pn, JM)) is the two quasi-particle creation (annihilation) operator coupled 
to the good angular momentum J with projection M, namely 



A\pn,JM)= Y. CI 



Jpmpjnirin prup'^nmn " 



(39) 



In the pn-RQRPA the commutator of two-bifermion operators is replaced by its ex- 
pectation value in the correlated QRPA ground state \Oq^pj^ > (renormalized quasiboson 
approximation). Therefore we have 

{OliPA\[Mpn,JM),A+{p'n',JM)]\0+p^)=6,,,6nn' x 

{l - — < O^p^|[a+ap]oo|0^p^ > - — < O^p^l [a+a„] oo|0^pa > \, 



Jl 



3k 



V 



pn,J'^ 



(40) 

with jp = J2jp + 1. Replacing \Orpa > in Eq. (^Of ) by the uncorrelated BCS ground state 
leads to the usual quasiboson approximation ( Vpn = 1), which violates the Pauli exclusion 
principle by neglecting the terms coming from the commutation rules of the quasi-particles. 
From Eqs. (^) one can derive the RQRPA equation 



1/2 A B \ 1/2 1/2 / X- 

\ -B -A I F" 



^7. 1)1/2 






(41) 



A,B 



jrm ym 



X^,Y^ 



The matrices A and B are given as follows: 

'A{pn, JM), \h, A^ip'n', JM) 



- 7"" 

r 






1/2 
pn,J^ \^RPA 



Q+ \ r)-l/2 

Ep + E.Jpp'6nn' - 2[ G{pn,p'n'; J){UpUnUp,Un' + UpUnUpmn') - 

F{pn,pn; J){UpVnUp'Vn' + VpUnVp'Un') ] DI'^j^DHIij^, 



T^J" 



B. 



pn,p'n' 



^pn,J^ V^RPA 



^RPA 



I ^p'n', 



A{pn, JM), [H, A{p'n\ JM) 

'^DI^u^J-DHIij^ [ G{pn,pn; J){UpUnVp>Vn> + VpVnUp'Un') 
F{pn,p'n] J){UpVnVp'Un' + VpUnUp'Vn') ]■ 



1/2 



(42) 



(43) 



Here, G{pn,p'n', J) and F{pn,p'n', J) are the particle-particle and particle-hole interaction 
matrix elements, respectively [Q. The coefficients Vpn^j^r are determined by solving numer- 
ically the system of equations |p|,ll| : 



V. 



pn,J'" 



1- 



2ip + 



E7-, 7/2 iv™ 1 2 

^pn',,/"^'^ l-"^ (pn',J"^')l 
-■- „' 



J"r', 



2jn + 



ET-, T/2|T7™ 1 2 

'^p'n,J'^'^ \'' (p'n,J"^')l • 



(44) 



v 
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The selfconsistent scheme of the calculation of the forward- (backward-) going free variational 
amplitudes X (V ), the excited energies f^™^ related to the ground state and the coefficients 
T^pn,.]-^ is a double iterative problem which requires the solution of coupled Eqs. (^) and 



Numerical treatment of the matrix elements (p0|)-(p3|) within the pn-RQRPA needs trans- 
formation of them to relative coordinates. After some tedious algebra one can obtain 



Kve= E (-y"+v+^+^(2:r + i)|f f i] 



X 



pnp' n 
J'^ irij^ m j- J' 



<p(l)^'(2);:r|/(ri2)r+r+(94,,(12)/(ri2)|n(l),n'(2);:r>x 

< O;^ II [c^,Cn'].j II J^nif >< rmf\rmi >< J'^TJii \\ [c^Cn\j \\ 0+ > . (45) 

In Eq. (^3]) Ofypg(12) represents the coordinate and spin dependent part of the nuclear two 
body transition operator for the 0i//3/3-decay and can be expressed in the form: 

Olypei'^2) = Hly^^_p{ru) + Hly^^_aAri2)cri2 + H^yp^_An2)Su. (46) 

The short-range correlations between the two interacting protons (p(l) and p'{2)) and neu- 
trons (^(1) and ^'(2)) are simulated by using the correlation function firu) in the non- 
antisymmetrized two-body matrix element (Eq. (P5|)). We adopt its following form: 

/(na) = l-e""^?2(l-6r^2) with a = 1.1 fm^ and 6 = 0.68 fml (47) 

The one-body transition densities entering Eq. (^) are given as follows: 

< .rruj II [c+Cn]j \\0i > _ (.) (.) ^, (.. ,.^^^ N./-r,W (ao\ 

/__ _ — [Up U^ v\(p„ j^) -\-Vp U^ r {pn,J^))y ^pnj^y \^°) 



V2J 



< Ojr II [c+Cn]j II J^'nif > 
V2J+1 



f„(/)^,(/)X"-f ^ I ^,(/)^,(/)y™/ ^,/x)(/) (Aq) 

[Up U„ ^ {pn,J^) ^ '^p ^n ^ {pn,J^))\J ^pn,J^- \^^) 



The index i (f) indicates that the quasiparticles and the excited states of the nucleus are 
defined with respect to the initial (final) nuclear ground state |0^ > (|0j >). As the two 
sets of intermediate nuclear states generated from the initial and final ground states are not 
identical we use the overlap factor 

^ 7+ I 7-+ ^ ~ ^T dT""^ 'Y"'f - v^^ V"^f \ V 

^ -'m^l-'m. ^ ~ Z^ y^ {pn,J^)^ {pn,J^) ^ {pn,J^)^ {pn,J^)J ^ 
pn 

(4'^<^ + ^i^^if^)i^M^ + ^i^^^i^^) (50) 

in definition (ESf) of the nuclear matrix elements. 
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IV. RESULTS AND DISCUSSION 

A. Calculation of the nuclear matrix elements 

The pn-RQRPA method has been apphed for the calculation of the SUSY Oz//?/5-decay 
nuclear matrix elements of the A = 48, 76, 100, 116, 128, 130, 136 and 150 nuclear systems. 
The following single particle model spaces have been considered in these cases: 
(i) For "^^Ca ^^ Ti decay the nuclear model comprises 13 levels: Os^ /j, Op^^m, Opg^, Is^^, 

0*^3/2' 0*^5/2' -'-Pl/2' -'^^3/2' 0/5/2' 0/7/2' ^S-^^g' ^97/2^ ^99/2' 

(ii) For "^^Ge — s>'^^ Se, ^^Se — >^^ Kr decays the model space comprises 12 levels: ls^,2, 0d^,2, 

0*^5/2' lPl/2' 1^3/2' 0/5/2' 0/7/2' 2^1/2' 1^3/2' 1*^5/2' Ofi'7/2, ^99/2' 

(iii) For ^^Zr -^^^ Mo, ^^^Mo ^^°° Ru and ^^^Cd ^^^^ Sn decays the model space com- 
prises 16 levels: Isy^, ^dy^, Oc^g/a, lPi/2, 1^3/2' O/5/2' O/7/2, 25^/2' 1^^3/2' 1^^5/2, 0^7/2' Oc/g/a, 

1/5/2' 1/7/2' 0/^9/2' O/iii/2- 

(iv) For ^2^Te ^^^^ Xe, ^^°Te ^^^° Xe and ^^^Xe ^^^^ Ru decays the model space com- 
prises 16 levels: 15^/2, 0c?3/2, Orfg/S' lPl/2' 1^3/2' O/5/2, O/7/2, 25^/2' 1^3/2' lc?5/2, 0^7/2' 0%2' 
2Pi/2' 2^3/2' 1/5/2' 1/7/2' 0/^9/2' O/in/2- 

(v) For ^^^Nd — i>i50 5*771 decays the model space comprises 20 levels: lSj^/2, 0^3/2, 0d^,2, 

lPl/2' 1^3/2' 0/5/2' 0/7/2, 2Si/2, lc?3/2, 1^^5/2' 05'7/2, ^99/2^ '^Pl/2^ '^P3/2^ I/5/2' I/7/2' ^K/2^ 
0^11/2' Oiii/2, 07^3/2. 

The single particle energies have been calculated with a Coulomb-corrected Wood-Saxon 
potential. For the two-body interaction we used the nuclear G-matrix calculated from the 
Bonn one-boson exchange potential. The single quasiparticle energies and occupation am- 
plitudes have been found by solving the BCS equations for protons and neutrons. Since 
our model spaces are finite the proton-proton and neutron-neutron pairing interactions have 
been renormalized according to Ref. |^S[. In the calculation of the pn-RQRPA equation 



we also renormalized the particle-particle and particle-hole channels of the G-matrix in- 
teraction by introducing two parameters Qp^ and f^^^, i.e. G{pn,p'n', J) -^ gppG{pn,p'n', J) 
and F{pn,p'n', J) -^ gphF{pn,p'n', J). The value adopted from our previous calculations 
|l0| , p!2|j36| , |39[| is gph = 0.8 and we discuss the stability of the nuclear matrix elements in 
respect to the gpp inside the expected physical range 0.8 < gpp < 1.2 . 

By the method outlined above we obtained the particular nuclear matrix elements of the 
two-nucleon mode Mgtn, -Mfn, M.gt' J^f', -Mt' and of the pion-exchange mode J^cT-kn, 
M-T-k-K for all above mentioned nuclei. Their values are listed in Tab. | for gpp = 1.0. It 
is worthwhile noticing that the two-nucleon matrix elements are considerably smaller in 
comparison with the pion-exchange ones. From that table one can see that the smallest 
nuclear matrix elements are associated with A=48 and 136 systems. We suppose that it 
is connected with the fact of the closed shell in ^^Ca and ^^^Xe. Sharp Fermi level usually 
offers weaker transitions due to the Pauli blocking effect. The largest nuclear elements are 
associated with A=150 and 100 systems. The difference between the Gamow- Teller matrix 
elements of A= 48 and A=150 systems is about an order of magnitude. We also stress that 
the one-pion exchange mode is disfavored in respect to the two-pion exchange mode not 
only because of a considerably smaller value of the structure coefficient {a^'^ = —0.044 ^ 
a^'^ = 0.20) but also due to the partial mutual cancellation of J^cr-in and A4t-1tt in Eq- 
(^) for all studied nuclear systems. By multiplying the particular nuclear matrix elements 
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with corresponding structure coefficient one obtain the full two-nucleon Ai'^^ , Aiy^ and 
pion-exchange Ai'"'^ matrix elements [see Eqs. (|30D-(|32D]. The A^?^ and Aiy^ have been 
calculated by using both the NR and the BM structure coefficients. By glancing at Tab. | we 
can see that the values of Ai"^^ are considerably larger comparing to Ai^^ and less sensitive 

to the chosen type of structure coefficients. Also the pion-exchange matrix elements A4'^^ 
dominate by a factor 5-7 over A^?-^ element for all studied nuclear systems. 

The sensitivity of the nuclear matrix elements A^?^, Ai^^ and Ai'^'^ to the details 
of the nuclear model and to the effect of short-range correlations are presented in Table 
^ The advantage of the pn-RQRPA in respect to the QRPA method is that there is no 
collapse of the pn-RQRPA solution within the physically acceptable region of the particle- 
particle interaction parameter Qpp (0.8 < Qpp < 1.2). We can see that the two-nucleon 
matrix element A4W is rather small and unstable in respect to the changes of Qpp within 
the discussed interval. Moreover, its value crosses zero for most nuclear systems. It means 
that if the deduced limits on A'^^^^^ associated with this nuclear matrix element only, the 
predictive power of the result would be rather small. The matrix elements Al|^ and Al'^^ 
are suppressed by the repulsion of the nucleon-nucleon interaction at short distances by the 
factor of about 5 and 3, respectively. Thus the large value of Al'^^ in comparison with Al^^ 
is predominantly due to the effect of the bosonization of tt^ ^ 7r+ + 2e^ vertex. We note 
that variations of the values of Al^^ and Al'^^ do not exceed 20 % - 30 % in respect to the 
their average values within the allowed gpp interval. Fig. |l| presents the dependence of the 
two-nucleon mode (Al?^ and Al|^) and the pion-exchange mode (Al'^^) calculated with 
and without short range correlations on the parameter Qpp for A=76 system. One can see 
the large suppression of the matrix elements due to the short range correlations effect as well 
as rather stable behavior of the results on Qpp. The inclusion of ground state correlations 
beyond the QRPA within the pn-RQRPA method stabilizes the behavior of the studied 
matrix elements as a function of Qpp and increases their predictive power, which is consistent 
with other studies [^ . 



B. Constraints on R-parity violation 

Here we shall analyze the constraints on R-parity violation using the current experimental 
lower half-life limits of the Oz//?/3-decay for the nuclei listed in Table |TT|. A subject of our 
interest is the dependence of the limits of three commonly used in literature quantities 
A'm, \'ni/{{mq/WOGeVf{mg/WOGeVy/^) and A'^/ (("^e/100 GeV)2(m^o/100 GeV)!/^) 
on the GUT scale scenarios of the MSSM. Moreover, we compare contributions to final 
limits coming from two-nucleon and pion-exchange modes and also study the relevance of 
the gluino and the neutralino exchange mechanisms. 

In the ffist step, we pay our attention to the constraints on free parameters of MSSM, 
particularly such SUSY parameters as tan/3 = — , mo, 'n^i/2, Aq and sign{fi). The allowed 
space for these parameters is determined by some natural restrictions. We use the condition 
that the mass eigenstates of SUSY particles cannot be imaginary and consider the dynamical 
electroweak symmetry breaking conditions (Eq. ([I9|) ). Further, we require that the i?7 
parameter associated with flavor changing neutral current processes fulfills relation (p2t). 
Using all the above assumptions, we have determined the excluded region for the parameters 
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tan/3 and rrio. It is shown in Figs, ^(a) and |](b) where corresponding symbols indicate three 
main sources of constraints. The caption EWSB indicates the points excluded by improper 
electroweak symmetry breaking, 6 — >■ s + 7 the points eliminated by FCNC constraints and 
v.e.v. - the points, for which some of the mass eigenstates become imaginary. One can 
observe that this region depends strongly on sgn(/i): for /i being positive, the allowed region 
of tan/3 and mo is considerably larger. We find e.g that for yU > and small values of tan/? 
(< 3) there is no restriction on rriQ while for fi < there exists a lower bound on mo of 
about 240 GeV. 

Further we have found that the constraints coming from the B -^ Xg'j decay are also 
especially sensitive to sgn(/i). In Figs, ^(a)-(b) the dependence of the R^ parameter on 
mo and mi/2 is shown, for both positive and negative values of /i. We note that the SM 
contribution to R7 is fixed for all the MSSM points at -0.188. The main SUSY contributions 
to this parameter come from charged Higgses and charginos and are shown on Figs. @(c)-(d) 
for negative fi. One can see instantly the large negatives contribution coming from charged 
Higgses which is responsible for exclusion of the corresponding low values of tuq and mi/2- 

We proceed with the determination of limits on 

A'lii, X[u/ {i^g/WO GeVfim-g/lOO GeVy/^) and A'^/ (("^e/100 GeV)2(m^o/100 GeV)!/^) 
from the non observation of double beta decay as a function of SUSY free parameters. For 
this purpose, we use the Oi^/5/5-decay half-life time ( p9D and the calculated nuclear matrix 
elements. We have found a very weak dependence of the quantities under discussion on the 
tan /3 and Aq SUSY parameters (tan /3 influences mainly the Higgs (and neutralino) sector, 
while Aq determines a very weak sfermion mixing). Therefore we present dependence on 
m,i/2 and mo for positive /z only. The results are drawn in Figs, ^(a)-(c) for all nuclei we are 
interested in. One can see that the currently strongest bound on these quantities is deduced 
from the A=76 system followed by A=128. The limit on the trilinear ^p breaking SUSY 
parameter X'^^-^^ becomes less stringent with the increasing value of mo (Fig. §(a)). Similar 
behavior is found for A'^/ ((m,-/100 GeV)2(mg/100 GeV)^/^^ also (Fig. |(b)). On the con- 
trary, the value of X'm/ f (mg/lOO GeV)^(mj^o/100GeV)^/^j decreases with the increase of 
m^o- The reason for such a behaviour is connected with the fact that the mass of selectron 
is growing faster than the mass of squark as the value of m,o is increased. In addition, the 
increase of its second power compensates for the modest increase of X[^i in respect to mo. 
This can be seen explicitely from the dependence of the of squark, selectron, gluino and the 
lightest neutralino masses on mo as drawn in Fig. ^. We note that gluino and neutralino 
are gauginos and therefore are insensitive to the common scalar mass at the GUT scale mo. 



In Tab. |Tl| we present the limits on ^p coupling constants X'm from the current lower 



bounds on the half-life time of the OuPP-decaj isotopes experimentally most promising. We 



consider two different scales for mo and m,i/2: 100 GeV and 1 TeV. By glancing at Tab. pi 
one finds out that the strongest limit on X'm comes from the ^^Ge isotope and is 5.25 x 10~^ 
and 1.8 x 10^^ for the considered values of mo and mi/2, respectively. 

The experimental constraints on the half-life of the Oi//3/3-decay are expected to be more 
stringent in future. It would require the recalculation of the corresponding limits on X'^^^. 
For the experimentalists convenience and in order to make our results more general we 
introduce the GUT constrained SUSY sensitivity parameter ^y^"^'^^^(mo, mi/2) of a given 
isotope Y, presented in Tab. |T|. This parameter incorporates the elements of both ^pSUSY 
and nuclear structure theory and is related with the limit on A'^^^^ as follows: 
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A'ni(mo,mi/2) < er^^(mo,mv2) x f o./3/e.p l • (51) 

In this way the reader is provided with an easy algorithm of predicting desired hmits with 
changing the experimental data. 

A detailed study of the GUT constrained SUSY scenario for the Oz//3/3-decay of '^^Ge 
offers the most stringent limit on R-parity breaking. At present this process is also the 
most perspective for the experimental detection in respect to the planned future experiment 
GENIUS |20]. In Figs. |^(a) and ^(b) we show the most stringent restriction on )k^^^ and 



'^111/ (("^e/100 GeV)^(m^o/100 GeV)^/^ j as a function of both mo and mi/2. In order to be 
more concrete we show the present bounds on \'^^^ together with those expected from the 
GENIUS Oi//3/3-decay experiment for a given set of mo and ?7ii/2 parameters in Tab. |^. One 
can see that the GENIUS experiment is expected to improve the limits on A'^^^^^ by about 
factor of 5 in respect to the current ones. Tab. |I^ contains also the calculated values of 
the corresponding masses of squark, selectron, gluino and lightest neutralino as well as the 
sensitivity parameter ^ncf^^ within the GUT constrained SUSY scenario. 

We also discuss the importance of different Oz//5/5-decay mechanisms on hadron and quark 
levels. The theoretical expression for half-life of -^pOz//5/5-decay in Eq. (|29D comprises the 
contributions from the two-nucleon mode and pion-exchange mode, which incorporate a 
different combination of ^p SUSY parameters tjt and rjps- In the previous subsection we 
have shown that the pion-exchange nuclear matrix elements dominate over the two-nucleon 
ones. However, it is not the necessary condition for the dominance of pion-exchange mode. 
In order to clarify this point we calculate the limit on )k^^^ by considering only one of these 
mechanisms at one time. The dependence of A'^^^^ on ttiq is drawn in Fig. ^ for ?7ii/2 = 100 
GeV and 500 GeV. We see that the pion-exchange mode offers a considerably stronger limit 
on \\^^ than the two-nucleon mode, which can be safely neglected. We note that the peak 
appearing in the curves presenting the two-nucleon mode is a consequence of interference 
between neutralino and gluino contributions to final amplitude. The curve representing the 
limit on A'^^^^^ from the pion-exchange mode is free of such instabilities. 

It remains to find out which of the neutralino and gluino ^pexchange mechanisms is the 
most important one for the Oi//3/9-decay process. In Fig. |^ we show the parameter \'^^^ as 
a function of mo for mi/2 equals to 100 GeV and 500 GeV by considering only one of the 
above ^pmechanisms. For ?7ii/2 = 100 GeV and ttlq larger than about 200 GeV the gluino 
mechanism is the dominant one. On the other hand, for mi/2 = 500 GeV the neutralino 
exchange mechanism becomes prevalent within the whole considered interval of ttiq. One 
concludes that none of these mechanisms is the most important in general. The answer to 
this problem depends on the details of the GUT constrained SUSY scenario and in particular 
on the values of mo and mi/2 parameters. 

V. CONCLUSIONS 

In conclusions, we have presented a detailed analysis of the R-parity violating trilinear 
terms contribution to the Oi//5/5-decay within the GUT constrained MSSM scenario. Both 
nuclear and particle physics aspects of this process have been discussed to some extent. 



We calculated relevant nuclear matrix elements within the realistic pn-RQRPA method. 
A comparison of nuclear matrix elements belonging to the two-nucleon mode and the pion- 
exchange mode has been performed for A = 48, 76, 82, 96, 100, 116, 128, 130, 136 and 150 
nuclear systems. We have found that the vr mechanisms are larger by factor 5-7 mostly due 
to the strong ^p vr^ ^ vr^ + 2e^ transition. We suppose it is the explanation why the pion- 
exchange mode Oi//3/3-decay is not favored for the mechanisms based on exchange of Majorana 
neutrinos having a different Lorentz structure and predicting weaker ^p tt" ^ vr"*" + 2e~ 
transition. Our studies show that ^pp ion-exchange mode nuclear matrix elements are less 
suppressed by the short-range correlation effects and are more stable in respect to the details 
of nuclear Hamiltonian than the two-nucleon matrix elements. 

In the MSSM part of calculations we limited the space of the free SUSY parameters using 
necessary conditions for proper electroweak symmetry breaking and inducing other limits 
based on current FCNC (Flavor Changing Neutral Currents) experiments. The exclusion 
plots for the SUSY parameters are presented and the sensitivity of the excluded region to 
the sign of /i parameter is manifested. A detailed study of the R-j parameter crucial for 
analysis of the FCNC B -^ Xs7 decay processes is also displayed and analyzed. We also 
present masses of squark, selectron, gluino and lightest neutralino in the GUT constrained 
MSSM, the knowledge of which is required for the analysis of the ]/lp parameters. 

Using the experimental lower bounds on the Oz//5/5-decay half-life we then deduced current 
constraints on the ^pMSSM parameters for different nuclei. The most restrictive constraints 
on the ^p Yukawa coupling constant X'^^-^ are found from the ^^Ge Oz//3/3-decay experiment 
performed by the Heidelberg- Moscow Collaboration [|l^ . For the common SUSY masses ttlq 
and mi/2 at the scale of 100 GeV and 1 TeV A'^^^^ they are: 

A'ni < 5.3 X 10"^ present, 1.1 x 10"^ (GENIUS) for mo = mi/2 = 100 GeV, (52) 
A'ni < 1.8 X 10"^ present, 3.7 x 10"^ (GENIUS) for tuq = mi/2 = 1 TeV. (53) 

The present limit on A'^^^^^ can be improved by a factor of 5 if the future Oz//3/3-decay experi- 
ment GENIUS is carried out and no signal about the Oi//3/3-decay is detected. 

Our studies have shown that the above limits are associated with the pion-exchange 
mode of the Oi//3/3-decay, which is the dominant mechanism for this process at the hadron 
level. Therefore, bearing in mind the above nuclear structure analysis we argue that the 
obtained limits do depend very insignificantly on the nuclear physics uncertainties. 

We have also dealt with the question which of the gluino and neutralino Oz//3/3-decay 
mechanisms is more important. We have shown that there is no unique answer to this 
problem and the dominance of any of them is bound with a different choice of the SUSY 
parameters rriQ and mi/2- 

Finally, we conclude that the Oi//3/3-decay imposes very restrictive limits on the important 
R-parity violating SUSY parameters also within the GUT constrained MSSM scenario. The 
obtained constraints on the trilinear ^pSUSY parameter A'^^^^^ are able to compete with those 
derived from near future accelerator experiments. 
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TABLES 



TABLE L Nuclear matrix elements (30 - ^3]) for the R-parity violating SUSY mode of the 
Ovfifi-decay for the experimentally most interesting isotopes calculated within the renormalized 
pn-QRPA with gpp = 1.0. AlxlO" implies that the matrix element should be divided by 10" to 
get the current numerical value. BM (NR) denotes that the nucleon structure coefficients of the 
bag model (non-relativistic quark model) have been considered. 

(/3/3)ojy — decay : 0^^ -^ 0^^ transition 

M. E. 48^„ 76q^ 82^^ 96zr ^^^Mo ^^<^Cd ^^^Te ^^"^Te ^^^Xe i^o^^ 

two-nucleon mode 

Mgt^n X 10^ 1.45 7.05 6.52 4.54 8.14 4.91 7.37 6.64 3.87 11.0 

Mf-n X 10^ -0.58 -2.48 -2.28 -1.67 -2.94 -1.78 -2.68 -2.43 -1.44 -4.09 

Mgt' X 10^ -0.21 -1.04 -0.99 -0.64 -1.17 -0.74 -1.05 -0.94 -0.54 -1.54 

Mf' X 10^ 0.92 3.76 3.53 2.46 4.41 2.79 4.00 3.61 2.14 6.05 

Mt' X 10^ -1.30 -2.38 -1.98 -2.68 -4.14 -1.88 -3.67 -3.46 -1.92 -6.34 

Ml^ (BM) -19.6 -95.9 -87.7 -63.7 -113. -65.4 -102. -92.0 -53.1 -154. 

Mf^ (BM) -0.32 4.71 5.01 1.02 3.15 3.10 2.78 2.19 1.29 2.55 

A^?^ (NR) -27.4 -129. -118. -85.1 -151. -89.3 -137. -123. -71.9 -207. 

Mf^ (NR) -1.06 2.22 2.80 -0.89 -0.07 1.28 -0.17 -0.52 -0.32 -2.13 

pion mode 

Mgt-Itt 0.25 1.30 1.23 0.77 1.43 0.92 1.25 1.10 0.61 1.85 

Mr-in -0.53 -1.02 -0.87 -1.11 -1.73 -0.78 -1.57 -1.48 -0.84 -2.70 

MGT-27r -0.26 -1.34 -1.26 -0.85 -1.52 -0.94 -1.40 -1.26 -0.74 -2.07 

MT-2n -0.31 -0.65 -0.57 -0.67 -1.05 -0.47 -0.99 -0.93 -0.54 -1.68 

7W^ -147. -625. -583. -428. -750. -436. -692. -627. -367. -1054. 
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TABLE II. Sensitivity of nuclear matrix elements A^|^, MY and M""^ for the A=48, 76, 
82, 96, 100,116,128,130,136 and 150 nuclear systems calculated within the renormalized pn-QRPA 
to of nucleon short-range correlations (s.r.c) and to the factor Qpp, renormalizing the parti- 
cle-particle interaction strength. In the calculations of nuclear matrix elements M^ and A4|^ 
within two-nucleon mode the nucleon structure coefficients of the non-relativistic quark model was 
adopted. 



nucleus 


Qpp 




Mf 


Mf 




M""^ 




no s.r.c. 


with s.r.c. 


no s.r.c. with 


s.r.c. 


no s.r.c. 


with s.r.c. 


^^Ca 


0.80 


-154. 


-31.2 


-22.1 


-0.85 


-440. 


-164. 




1.00 


-138. 


-27.4 


-20.3 


-1.06 


-392. 


-147. 




1.20 


-119. 


-23.0 


-18.1 


-1.27 


-335. 


-126. 


76Ge 


0.80 


-711. 


-143. 


-90.8 


3.67 


-2020. 


-686. 




1.00 


-646. 


-128. 


-80.5 


2.22 


-1831. 


-625. 




1.20 


-581. 


-113. 


-78.2 


0.81 


-1645. 


-564. 


82^6 


0.80 


-638. 


-131. 


-79.6 


4.13 


-1832. 


-638. 




1.00 


-581. 


-119. 


-74.2 


2.80 


-1667. 


-583. 




1.20 


-525. 


-105. 


-68.8 


1.52 


-1503. 


-529. 


82Zr 


0.80 


-490. 


-97.0 


-66.3 


0.25 


-1383. 


-479. 




1.00 


-437. 


-85.1 


-61.1 


-0.89 


-1228. 


-428. 




1.20 


-358. 


-73.4 


-55.9 


-1.85 


-1077. 


-376. 


iooMo 


0.80 


-850. 


-170. 


-112. 


1.86 


-2408. 


-832. 




1.00 


-764. 


-151. 


-104. 


-0.07 


-2155. 


-750. 




1.20 


-678. 


-131. 


-95.5 


-1.80 


-1908. 


-668. 


"6cd 


0.80 


-480. 


-98.2 


-61.5 


2.05 


-1365. 


-474. 




1.00 


-440. 


-89.3 


-57.5 


1.28 


-1250. 


-436. 




1.20 


-401. 


-80.3 


-53.5 


0.55 


-1135. 


-398. 


I28re 


0.80 


-778. 


-155. 


-103. 


1.81 


-2212. 


-767. 




1.00 


-696. 


-137. 


-95.1 


-0.17 


-1976. 


-692. 




1.20 


-615. 


-118. 


-87.5 


-2.08 


-1743. 


-615. 


laoj^g 


0.80 


-706. 


-141. 


-93.6 


1.37 


-2008. 


-697. 




1.00 


-631. 


-123. 


-86.7 


-0.52 


-1790. 


-627. 




1.20 


-556. 


-106. 


-79.9 


-2.36 


-1573. 


-556. 


136^6 


0.80 


-418. 


-83.5 


-55.0 


1.04 


-1191. 


-413. 




1.00 


-369. 


-71.9 


-50.8 


-0.32 


-1048. 


-367. 




1.20 


-318. 


-59.7 


-46.6 


-1.74 


-901. 


-317. 
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150 



Nd 



0.80 


-1185. 


-234. 


-160. 


0.52 


-3360. 


-1167 


1.00 


-1066. 


-207. 


-149. 


-2.13 


-3013. 


-1054 


1.20 


-949. 


-180. 


-138. 


-4.61 


-2677. 


-943 



TABLE III. Upper limits on the lepton number non-conserving parameter X'm deduced from 
the experimental lower limits of the Ou (3 f3-decay half-life time T^J^ ^'^^^(Y) for the nuclei studied 
in this work. The MSSM SUSY parameters ttiq and m.1/2 are limited to two cases: 100 GeV and 1 
TeV. According to Eq. (51) £,y^^^^ (''^011^1/2) denotes the sensitivity of a given nucleus Y to the 



A'^]^^ parameter. Gqi is the integrated kinematical factor for 0^^ -^ 0^^ 


transition. 




mo = m.1/2 = 100 GeV 
nucleus Gqi x lO^^y rf/f -^^^(Y) [y] Cf^^*^ A^ii 


mo = mi/2 = 

eMSSM 


= 1 TeV 

■^111 



^^Ca 


803. 


'^Ge 


7.93 


825e 


35.2 


82Zr 


73.6 


iooMo 


57.3 


ii^Cd 


62.3 


i28Te 


2.21 


130^g 


55.4 


i36Xe 


59.1 


150 JY^ 


269. 



9.5 X 10^1 ^ 

1.1 X 10^5 I 
2.7 X 10^2 I 

3.9 X 10^9 I 

5.2 X 10^2 I 
2.9 X 10^2 I 
7.7 X 1024 I 
8.2 X 10^1 I 
4.2 X 10^3 I 
1.2 X 10^1 I 



5.90 X 
9.56 X 
6.84 X 
6.27 X 
5.26 X 
6.83 X 
1.23 X 
5.79 X 
7.44 X 
2.95 X 



10-4 
10-4 
10-4 
10-4 
10-4 
10-4 
10-3 
10-4 

10-4 
10-4 



1.9 X 10-3 
5.2 X 10-4 
1.7 X 10-3 
7.9 X 10-3 

1.1 X 10-3 
1.6 X 10-3 
7.4 X 10-4 
1.9 X 10-3 

9.2 X 10-4 
1.6 X 10-3 



0.202 
0.327 
0.234 
0.215 
0.180 
0.234 
0.422 
0.198 
0.254 
0.101 



0.65 
0.18 
0.58 
0.27 
0.38 
0.57 
0.25 
0.66 
0.32 
0.54 
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TABLE IV. The supersymmetric parameters with corresponding hmits on X'^^^ derived from 
the best presently available experimental limit on half-life of ^^Ge Oz//3/3-decay: T^f^C^^Ge) 
> 1.1 X 10^^ y [19]. Consequences of the expected half- life limit to be reached in the GENIUS 
experiment {T^^^^^^(^^Ge) > 6 x 10^''' [18]) are also shown. The tuq, mi/2 tan(/3), ^o and sgn(/x) 
are the MSSM parameters, rriq, rrie, rrig and mi are the masses of squark, selectron, gluino and 
lightest neutralino, respectively. The £,^q^^ parameter is defined by Eq. (51). 









tan(/?) = 


3, Ao = 500 GeV, sign{fi) = +1 






mo 


"1-1/2 


nig 


rrie 


rUg 


m^,i 

Xo 


tMSSM 


^111 


^111 


[GeV] 


[GeV] 


[GeV] 


[GeV] 


[GeV] 


[GeV] 




present 


GENIUS 


100 


100 


250.9 


107.2 


261.9 


25.8 


9.56 


5.25 X lO"'' 


1.09 X lO-'' 


100 


500 


1136.0 


217.7 


1290.0 


208.0 


1.04 


5.73 X 10-3 


1.19 X 10-3 


100 


1000 


2281.0 


399.5 


2595.0 


420.2 


4.99 


2.74 X 10^2 


5.67 X 10-3 


500 


100 


548.1 


501.5 


298.7 


36.8 


1.15 


6.31 X 10^3 


1.31 X 10-3 


500 


500 


1237.0 


536.1 


1290.0 


208.3 


5.76 


3.16 X 10^2 


6.54 X 10-3 


500 


1000 


2333.0 


632.2 


2595.0 


420.2 


1.23 


6.76 X 10-2 


1.40 X 10-2 


1000 


100 


1025.0 


1001.0 


317.7 


39.8 


4.25 


2.33 X 10-2 


4.83 X 10-3 


1000 


500 


1511.0 


1019.0 


1407.0 


208.9 


1.57 


8.62 X 10-2 


1.78 X 10-2 


1000 


1000 


2490.0 


1072.0 


2596.0 


420.3 


3.27 


1.80 X 10-1 


3.72 X 10-2 
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FIGURES 

FIG. 1. The nuclear matrix elements of the two-nucleon mode (A^? , A^i ) and 

pion-exchange mode {Mf^^) of the ^pSUSY contribution to the Oi^/3/3-decay of '^^Ge calculated 
within the pn-RQRPA. They are functions of the particle-particle interaction strength gpp with (a) 
and without (b) the two-nucleon short-range coorelations. 

FIG. 2. The constraints imposed by the B -^ Xgj process, dynamical electroweak symme- 
try breaking condition (ESWB) and vacuum expectation values condition on the tan/3 and uiq 
SUSY parameters. The excluded region is indicated by the corresponding symbols. We assume 
TTT-x /2 = 200 GeV and Aq = 500 GeV. Two cases of n negative (a) and fi positive (b) are shown. 
The allowed space of parameters is considerably larger for positive //. 

FIG. 3. The -R7 parameter determining the MSSM contributions to the flavor changing neu- 
tral current B -^ Xg'j decays is drawn against mo and m^^. The coefficient was calculated for 
tan/3 = 3, Aq = 500 and fi either positive (a) or negative (b). The charged Higgs (c) and chargino 
(d) contributions to Rj for negative fi are also presented. 

FIG. 4. The limits on (a) X[,, (b) (^,/iooGcV)^(m,/iooGeV)i/^ 

and (c) 7 — ,-,^„^ ^.^9/" — Mr^r^r^ ^T^^/9 dcduccd from the experimental lower bound on the half-life 

^ ' (mg/100GcV)^(m-o/100GeV)^''^ -^ 



X" 



of the Oz^/3/3-decay of different nuclei are plotted as a function of ruQ. Other free parameters are 
fixed as follows: Aq = 500 GeV, mi/2 = 500 GeV, tan/3 = 3 and fi > (see text for details). The 
needed nuclear matrix elements have been calculated within the pn-RQRPA. 

FIG. 5. The masses of super symmetric particles: rriq (squark), rUg (gluino), mg (selectron) 
and m^ (the lightest neutralino) as a function of niQ within the GUT constrained minimal super- 
symmetric standard model with R-parity breaking. Other free parameters are fixed: mi/2 = 100 
GeV, Aq = 500 GeV, tan/3 = 3 and /x > 0. 



FIG. 6. The deduced limits on (a) Aiii and (b) -, — ,-,^„^ ,r^9/"^ — ;,^^^ ■,rM/9 from the exper- 

^ ' J--'--'- ^ ' (mg/100 GeV)^(m-o/100GcV)-^'^ 

imental lower bound on the half-life of the Oz^/3/3-decay in '^^Ge plotted as a function of ttiq and 
"T'i/2- Other parameters are fixed as in the previous figures. 

FIG. 7. The deduced from the lower bound on the half-life of the Oz//3/3-decay in ^^Ge limit on 
X'lii plotted as a function of mo- Two-nucleon mode (2n.-mode) and pion-exchange mode (vr-mode) 
are considered separately. The parameters mi/2 have been chosen to be 100 GeV or 500 GeV. Other 
parameters are fixed as in the previous figures. The dominance of the pion-exchange mode is easy 
to be noticed. 
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FIG. 8. The limit on A'^^^^ deduced from the experimental lower bound on the half-life of the 
Oz^/?/3-decay of ^^Ge against uiq. The gluino (g-mech.) and neutralino (x-niech.) mechanisms in 
the R-parity violating scenarios are considered separately for two values of mi/2 (100 GeV and 
500 GeV). Other parameters are the same as in Fig. 5. One can observe that for m^rt = 100 
GeV and a small value of rriQ up to 200 GeV the neutralino contribution dominates (gives stronger 
constraints), while for larger values of rriQ gluino mechanism becomes more important. 
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